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ABSTRACT 

We present submillimeter observations for 136 of the 370 X-ray sources detected in the 1 Ms exposure of 
the Chandra Deep Field North. Ten of the X-ray sources are significantly detected in the submillimeter. 
The average X-ray source in the sample has a significant 850 jjlth flux of 1.69 ± 0.27 mJy. This value 
shows little dependence on the 2 — 8 keV flux from 5 x 10^^^ erg cm^^ s^^ to lO^^'' erg cm^^ s^^. The 
ensemble of X-ray sources contribute about 10% of the extragalactic background light at 850 /im. The 
submillimeter excess is found to be strongest in the optically faint X-ray sources that arc also seen at 
20 cm, which is consistent with these X-ray sources being obscured and at high redshift {z > 1). 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: 
active 



and mostly arose from the z > 1.5 and spectroscopically 
unidentified sources. 

In this paper we present submillimeter observations of 
the X-ray sources detected in the 1 Ms exposure of the 
Chandra Deep Field North (CDF-N). These data cover 
more than three times the area of the Barger et al. (2001a) 
submillimeter/X-ray data. Furthermore, the 1 Ms hard X- 
ray data go about a factor of 10 deeper than the 100 ks 
SSA13 data, providing a much larger sample with a wide 
dynamic range in X-ray flux. The 1 Ms X-ray source cat- 
alog presented in Brandt et al. (2001) contains 370 signif- 
icantly detected point sources. Source positions are accu- 
rate to within k, 0.6 — 1.7", depending on the off-axis angle 
and the number of detected counts. For our analysis, we 
adopt from Table 3 of Brandt et al. (2001) the effective 
photon index (F) values, which are based on the ratio of 
the counts in the hard and soft bands, and the observed 
frame hard band flux values. The R magnitudes and red- 
shifts are taken from the 1 Ms X-ray follow-up catalog of 
Barger et al. (2001b), which gives redshift identifications 
for 66 of the 136 X-ray sources with submillimeter obser- 
vations. 

2. SUBMILLIMETER OBSERVATIONS 

SCUBA jiggle map observations at 850 /im were taken 
during observing runs in 2000 April and 2001 March and 
May. The maps were dithered to prevent any regions of 
the sky from repeatedly falling on bad bolometers. The 
chop throw was fixed at a position angle of 90 deg so 
that the negative beams would appear 45" on either side 
east-west of the positive beam. Regular "skydips" (Light- 
foot et al. 1998) were obtained to measure the zenith at- 
mospheric opacities, and the 225 GHz sky opacity was 
monitored at all times to check for sky stability. Point- 
ing checks were performed every hour during the obser- 
vations on the blazars 0923-^392, 0954-f685, 1144+402, 
1216-F487, 1219-^285, 1308+326, or 1418+546. The data 
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1. INTRODUCTION 

Most of the energy density of the extragalactic X-ray 
background (XRB) resides in the hard X-ray band above 
2 keV. Early 100 — 300 ks Chandra observations, in com- 
bination with ASCA observations at the bright flux levels, 
resolved > 60% of the 2 — 10 keV XRB into discrete sources 
(Mushotzky et al. 2001; Giacconi et al. 2001; Garmire et 
al. 2001a; Tozzi et al. 2001). Now two unprecedented 1 Ms 
Chandra exposures (Brandt et al. 2001; Rosati et al. 2001) 
have essentially completely resolved the remainder of the 
XRB at these hard energies. 

Many of the hard X-ray sources in the ultradeep Chan- 
dra data appear to be highly absorbed systems with large 
hard to soft X-ray flux ratios. In these sources the rest- 
frame soft X-ray through near-infrared radiation is re- 
processed by dust and gas surrounding the central active 
galactic nucleus (AGN) , and the reradiated energy appears 
in the far- infrared (FIR). At high redshifts (z > 1) this 
FIR radiation is redshifted to the submillimeter and can 
be directly measured with the SCUBA camera (Holland 
et al. 1999) on the 15 m James Clerk Maxwell Telescope. 

Early searches for submillimeter counterparts to Chan- 
dra X-ray sources yielded mixed results because of the 
small sample sizes: Fabian et al. (2000; A2390 and A1835) 
and Hornschemeier et al. (2000, 2001; Hubble Deep Field 
North and its vicinity) found no submillimeter counter- 
parts to their hard X-ray sources, while Bautz et al. (2000; 
A370) found submillimeter counterparts to both of theirs. 
Barger et al. (2001a) used submillimeter data that cov- 
ered the entire 57 arcmin^ SSA13 Chandra X-ray field of 
Mushotzky et al. (2000) to look for overlaps between the 
hard X-ray and submillimeter populations. Although only 
one hard X-ray source was significantly detected in their 
submillimeter data, these authors found that the error- 
weighted sum of the submillimeter fluxes of all 20 hard X- 
ray sources in their sample was signiflcant at 20 ± 6 mJy 



were calibrated using jiggle maps of the primary cali- 
bration source Mars or the secondary calibration sources 
CRL618, OH231.8-I-4.2, or 16293-2422. SubmiUimeter 
fluxes were measured using beam- weighted extraction rou- 
tines that include both the positive and negative portions 
of the beam profile. 

The data were reduced in a standard and consistent 
way using the dedicated SCUBA User Reduction Facility 
(SURF; Jcnness & Lightfoot 1998). The new data were 
combined in the reduction process with the jiggle maps 
previously obtained by Barger, Cowie, & Richards (2000). 
Due to the variation in the density of bolometer samples 
across the maps, there is a rapid increase in the noise levels 
at the very edges, so the low exposure edges were clipped. 

The SURF reduction routines arbitrarily normalize all 
the data maps in a reduction sequence to the central pixel 
of the first map; thus, the noise levels in a combined im- 
age are determined relative to the quality of the central 
pixel in the first map. In order to determine the absolute 
noise levels of our maps, we first eliminated the ^ 3cr real 
sources in each field by subtracting an appropriately nor- 
malized version of the beam profile. We then iteratively 
adjusted the noise normalization until the dispersion of the 
signal-to-noise ratio measured at random positions became 
~ 1. The noise estimate includes both fainter sources and 
correlated noise. 

3. SUBMILLIMETER FLUX MEASUREMENTS 

The SubmiUimeter fluxes of the X-ray sources were mea- 
sured in an automated fashion. We first measured the 
SubmiUimeter fluxes at the positions of the 20 cm (> 3cr) 
sources from Richards (2000) that fell on our SCUBA 
maps. Any radio source that we detected above the 3cr 
level at 850 ^m was included in a "detection list". We 
next measured all > 5cr SubmiUimeter sources that did 
not have a radio counterpart; we identified two and added 
them to the list. Finally, we added all the Chapman et al. 
(2001) SubmiUimeter measurements (made in photometry 
mode) of radio sources that did not lie within our SCUBA 
jiggle maps. 

After compiling the detection list, we compared the co- 
ordinates of the X-ray sources with the coordinates of the 
sources in the detection list. If an X-ray source was found 
to lie within 3" of a source in the detection list, we identi- 
fied the list source as the counterpart to the X-ray source. 
All of the radio-|-submillimeter source counterparts had co- 
ordinates within 1.5" of the X-ray source positions. Once 
the cross-identifications were made, the sources in the de- 
tection list were removed from the SCUBA maps since the 
wide and complex beam patterns of the bright SCUBA 
sources can produce spurious detections at other positions. 
SubmiUimeter fluxes and uncertainties were then measured 
at all the unassigned X-ray positions using a recursive loop 
which removed X-ray sources that were detected in the 
SubmiUimeter above the 3cr level in a descending order 
of signiflcance prior to remeasuring the fluxes. This pro- 
cedure again avoids multiple or spurious detections of a 
single bright SubmiUimeter source. 

Of the 370 Chandra X-ray sources in the CDF-N sample, 
136 have 850 yum flux measurements with uncertainties less 
than 5 mJy and 109 have uncertainties less than 2.5 mJy. 
Five flux measurements were drawn from the Chapman et 
al. (2001) sample, including a signiflcant 15.7 ± 2.4 mJy 



detection. While we did not reproduce the SubmiUime- 
ter detection of one of the Chapman et al. sources which 
feU within our jiggle maps (VLA J123624-H621743 was ob- 
served by Chapman et al. in poor weather conditions and 
found to be a marginal 3cr detection (see their Table 1); 
the source was also detected at the 3.7tT level in the scan- 
map of Borys et al. 2001), the very high significance of the 
15.7 mJy source suggests that it is real. Complete exclu- 
sion of the Chapman et al. data has no significant effect 
on our subsequent analysis. 

In total, we significantly (> Scr) detect in the submil- 
limeter 10 of the 136 X-ray sources. The average submil- 
limeter flux per X-ray source is 1.30 mJy, and the error- 
weighted average is 1.69 ±0.16 mJy. The flnal SubmiUime- 
ter measurements are presented in the CDF-N 1 Ms X-ray 
follow-up catalog of Barger et al. (2001b) and are plotted 
versus hard X-ray flux (squares) in Figure 0. We include 
in the figure the SSA13 hard X-ray sources (circles) from 
Barger et al. (2001a) to better sample the bright, hard 
X-ray flux end. 
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Fig. 1. — SubmiUimeter versus hard X-ray (2 — 8 keV) fluxes 
for the X-ray sources in the CDF-N (squares) and SSA13 (circles) 
fields with SubmiUimeter measurements. Sources with > 20/iJy 
radio fluxes are denoted by filled symbols, and sources with > 
Scr SubmiUimeter fluxes are denoted by a second large symbol. 
Thirty of the 46 sources to the left of the dashed vertical line at 
4 X 10~^® erg cm~'^ s~'^ (including the two significant SubmiUime- 
ter sources that are also radio sources) are plotted at their hard 
X-ray flux limits. Two of these have hard X-ray flux limits less 
than 10~^^ erg cm~^ s~^ and are shown at a nominal value of 
1.1 X 10~^® erg cm~^ s~^. Only one source to the right of the verti- 
cal line has only a hard X-ray flux limit (1.25 X 10"-*^^ erg cm~^ s~^)- 

In order to test the significance of our results, we re- 
peated our SubmiUimeter measurement procedure for a 
large number of simulations. In the simulations we offset 
the X-ray sample a random amount in the RA and Dec 
directions and then measured the SubmiUimeter fluxes at 
these random positions. The average flux per source at 
these random positions was found to be accurately zero, 
but the dispersion from field to field was found to be 
0.27 mJy rather than the nominal 0.16 mJy, suggesting 
the presence of some systematic error. We use this larger 
uncertainty subsequently. The 95 percent confidence level 
is 0.45 mJy, and the 98 percent confidence level is 0.6 mJy, 
leaving little doubt that there is a significant excess 850 /im 
flux associated with the X-ray sample. Similarly, the ran- 
dom samples have an average of 1.2 3ct sources per field 
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and 95 and 98 percent confidence levels of 4 and 5 sources, 
respectively, so most of the significantly detected submil- 
limeter sources are real. 

4. CONTRIBUTION TO THE SUBMILLIMETER EEL 

The resolution of the FIR/submillimeter extragalactic 
background light (EBL) by the DIRBE and FIRAS ex- 
periments on COBE revealed that the FIR/submillimeter 
EBL has approximately the same integrated energy den- 
sity as the optical EBL. Thus, the absorption and reradia- 
tion of light by dust in the history of galaxy formation and 
evolution is extremely important; however, in order to re- 
construct the star formation and accretion histories of the 
Universe, we need to be able to differentiate between the 
FIR/submillimeter contributions made by starbursts and 
those made by AGN. 

About 20 — 30 percent of the 850 /^m EBL has been 
resolved into discrete sources by blank field SCUBA sur- 
veys (e.g., Barger, Cowie, & Sanders 1999). Blain et al. 
(1999) estimated from a SCUBA survey of massive clus- 
ter lenses which amplify distant submillimeter sources that 
more than 80 percent of the 850 fim EBL is now resolved. 
The ultraluminous SCUBA sources are believed to be the 
high redshift analogs of the local ultraluminous infrared 
galaxies (ULIGs; Sanders & Mirabel 1996), but there is 
an ongoing debate even as to whether the local ULIGs 
are dominantly powered by star formation or by AGN 
activity. Modelling suggests that ^ 10 — 20 percent of 
the FIR/submillimeter EBL may be contributed by galax- 
ies containing bright AGN (Almaini, Lawrence, & Boyle 
1999; Gunn & Shanks 2001). It is difficult to spectroscop- 
ically look for AGN signatures in the SCUBA sources — 
though in a few cases such signatures have been detected 
(e.g., Ivison et al. 1998; Barger et al. 1999) — due to the 
extremely faint nature of the vast majority of the opti- 
cal counterparts. However, since most of the flux density 
from the Chandra X-ray sample is expected to arise from 
AGN activity, we may use the present data to place an 
upper limit on the AGN contribution to the 850 /im EBL. 
(A submillimeter source detected in X-rays is not neces- 
sarily solely powered by accretion onto the AGN, as there 
may also be dust obscured starburst activity in the galaxy 
which contributes to the submillimeter flux.) 

Barger et al. (2001a) multiplied the hard XRB mea- 
surement from Vecchi et al. (1999) by the ratio of the total 
850 fj,m. flux to the total hard X-ray flux of the SSA13 hard 
X-ray sample to estimate the hard X-ray contribution to 
the 850 /im EBL. They determined the contribution to be 
about 10 percent. However, their approach assumed that 
a fixed submillimeter to X-ray flux ratio applies to all of 
the X-ray sources. We can quantitatively see that this is 
not a good assumption from Table ||, where we list the 
average X-ray flux and the average submillimeter flux and 
uncertainty for the CDF-N-I-SSA13 X-ray sources in each 
of four hard X-ray flux bins. Since the average submil- 
limeter flux is consistent with being constant as a function 
of X-ray flux, the ratio of the submillimeter to X-ray flux 
rises with decreasing X-ray flux. 

A better approach may be to translate the measured 
average 850 fj,m. flux per X-ray source to an X-ray source 
contribution to the 850 ^va EBL. The SCUBA jiggle maps 
cover an area of 194 arcmin^. Excluding the Chapman et 



al. (2001) sources, which lie outside this area, there are 
131 X-ray selected sources with an average flux per source 
of 1.27 ± 0.27 mJy, where the uncertainty is again based 
on the random samples described in § 0. This translates 
to an 850 ^^m EBL of 3.1 ± 0.66 x 10^ mJy deg-^, which 
is 10 ± 2 percent of the total 850 /im EBL, if we adopt the 
3.1 X lO'' mJy deg'^ value of Puget et al. (1996), or 7 ± 2 
percent of the total, if we adopt the 4.4 x 10^ mJy deg~^ 
value of Fixsen et al. (1998). Since our results indicate 
that the integrated 850 /xm EBL is proportional to the 
total number of X-ray sources, the overall contribution 
to the EBL may increase with additional weaker X-ray 
sources. However, the rapid convergence of the X-ray 
number counts below 10~^^ erg cm~^ s~^ (Garmire et al. 
2001b) suggests that this correction would not be large. 



5. SOURCE PROPERTIES 

If the submillimeter flux represents the reprocessing of 
the radiated energy in the more obscured AGN, then we 
may expect to see a dependence of the submillimeter ex- 
cess on X-ray hardness. However, an opposing tendency 
is introduced by redshift: 850 fim source detections are 
primarily expected to be at 2; > 1 because of the strong 
negative if-corrections at this wavelength that nearly com- 
pensate for cosmological dimming, while high redshift X- 
ray sources will appear softer in the observed frame than 
in the rest frame. 
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Fig. 2. — Submillimeter flux versus F for the X-ray sources in the 
CDF-N (squares) and SSA13 (circles) fields with hard X-ray fluxes 
above 5.0 X 10~^^ erg cm~^ s~^. Sources with > 3a radio fluxes 
are denoted by filled symbols, and sources with > 3cr submillimeter 
fluxes are denoted by a second large symbol. 

In Fig. g we plot 850 fiui flux versus effective photon 
index for the nearly complete CDF-N hard X-ray sample 
(squares) above 5.0 x 10~^^ erg cm~^ s~^ (2 — 8 keV), 
where the photon indices are well deflned. We include 
on the flgure the SSA13 hard X-ray sources (circles) from 
Barger et al. (2001a) to better sample the bright X-ray 
flux end. Of the eight significantly detected submillimeter 
sources which are present in this more restricted sample, 
six (including the SSA13 source) are clustered at low F 
values, which suggests they are obscured systems, while 
the remaining two (including a non-radio CDF-N source) 
have quite high F values, which suggests they are unob- 
scured systems or obscured systems at high redshift. It 
should be kept in mind that we expect there to be one 



or two spurious detections in the sample, so these objects 
could be false positives. 

A more robust result can be obtained by consider- 
ing the average properties of the sources. Dividing the 
CDF-N-I-SSA13 combined sample into hard sources with 
r < 1 and soft sources with F > 1, we find that 
the hard sources have an average submillimeter flux of 
1.77 ± 0.21 mJy, which is substantially larger than the 
value of 0.89 ± 0.24 mJy in the soft sample. This suggests 
that most of the submillimeter flux is indeed coming from 
the obscured sources, although even the soft sources have 
a significant submillimeter excess. 

In Fig. we plot redshift versus submillimeter flux for 
the unrestricted CDF-N (squares) and SSA13 (circles) 
samples. Sources with spectroscopic redshifts are sepa- 
rated from sources without by two dotted lines at z = 
and at z = 3.5. Below the z = line, we plot the spectro- 
scopically unidentified sources with magnitudes R < 24.1, 
which corresponds to the / < 23.5 definition of Barger et 
al. (2001a) for optically bright galaxies. These sources, al- 
though bright enough for spectroscopic identification, have 
not yet been observed; they will likely populate the same 
region of the plot as the low redshift sources. Above the 
z = 3.5 line, we plot the spectroscopically unidentified 
sources with magnitudes R > 24.1. These sources are 
generally too optically faint for spectroscopic identifica- 
tion. All but two of the significantly detected submillime- 
ter sources fall into this region of the plot, which is consis- 
tent with the conclusion of previous SCUBA surveys that 
most submillimeter sources are optically faint. The two 
significantly detected submillimeter sources which are not 
optically faint have been observed spectroscopically: the 
spectrum for the z — 1.013 source, which is also a ra- 
dio source, shows strong Balmer absorption features and 
[Ne III] in emission, so it appears to be a reliable identi- 
fication, whereas the spectrum for the z = 0.555 source, 
which is not also a radio source, looks relatively normal, 
and thus its identification with the submillimeter source 
could be spurious. 

Again it is instructive to consider the average submil- 
limeter properties of the X-ray sources, this time accord- 
ing to redshift bin. The average submillimeter flux is sig- 
nificant in both the z = to 1 (0.80 ± 0.21 mJy) and 
the z = 1 to 3.5 (1.09 ± 0.35 mJy) spectroscopic red- 
shift bins. Moreover, the average submillimeter flux of the 
unidentified optically faint X-ray sources is highly signifi- 
cant (1.68 ± 0.19 mJy). The average submillimeter flux of 
the unidentified optically bright X-ray sources is not signif- 
icant (0.72 ±0.56 mJy) but is consistent with these sources 
populating the same region of the plot as the sources with 
spectroscopic redshifts. 

Figure ^ shows how the unidentified optically faint X- 
ray sources tend to separate according to 850 /im flux into 
sources with radio counterparts and sources without. In 
fact, all but one of the significant SCUBA sources with- 
out a spectroscopic identification have radio counterparts. 
This result is again consistent with previous analyses of 
SCUBA data which foimd that targeted SCUBA surveys of 
optically faint radio sources are a far more efficient means 
of identifying SCUBA sources than random blank field sur- 



veys (Barger et al. 2000; Chapman et al. 2001). 

Quantitatively, the average submillimeter flux of the 
optically faint X-ray sources with radio counterparts is 
2.86 ± 0.28 mJy, which is significant at above the lOcr 
level. In contrast, the average submillimeter flux of the 
optically faint X-ray sources without radio counterparts is 
not significant (0.74 ± 0.25 mJy). 

Figure H illustrates one additional property we can learn 
about the X-ray sources with significant submillimeter de- 
tections: their millimetric redshifts (Carilli & Yun 2000; 
Barger et al. 2000). We have overlaid on the data three 
curves which show the submillimeter flux dependence of a 
redshifted 20 /iJy, 40 /iJy, or 100 /iJy radio source whose 
spectral energy distribution is assumed to be like that of 
the prototypical ULIG Arp 220. From these curves we 
can see that the unidentified optically faint sources that 
are significantly detected in the submillimeter are likely to 
fall in the redshift range between z = 1 and 3. 
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Fig. 3. — Redshift versus submillimeter flux for the X-ray sources 
in the CDF-N (squares) and SSA13 (circles) fields. Sources with 

> 3(T radio fluxes are denoted by filled symbols, and sources with 

> 3(T submillimeter fluxes are denoted by a second large symbol. 
Sources without spectroscopic identifications are plotted at either 
the bottom of the figure (the optical counterparts have magnitudes 
R < 24.1) or at the top of the figure (the optical counterparts have 
magnitudes R > 24.1), separated from the spectroscopically identi- 
fied data by the dotted lines. Overlaid on the data arc three dashed 
curves which show the expected tracks of a redshifted 20 fj,3j (left- 
most curve), 40 fJ.Jj, and 100 /.tjy radio source, assuming a spectral 
energy distribution like that of the ultraluminous infrared galaxy 
Arp 220. 
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Table 1 
Average Fluxes for the CDF-N+SSA13 Sample in Four X-ray Flux Bins 

/(2 - 8 keV) Bin Average /(2 - 8 keV) Average /(850 fim) 

(10^^® ergs cm"2 s"^) (10"^'^ ergs cm"^ s"^) (mjy) 



5-10 

10-30 

30 - 100 

100 - 500 



7.5 
20 

44 
180 



1.7 ±0.32 

1.8 ±0.29 
1.3 ±0.27 

0.25 ± 0.42 
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